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"^s^" ■ We consider N bosons occupying a discrete set of single-particle quantum states in an isolated 

' trap. Usually, for a given excitation energy, there are many combinations of exciting different number 

, of particles from the ground state, resulting in a fluctuation of the ground state population. As a 

^s^j . counter example, we take the quantum spectrum to be logarithms of the prime number sequence, 

and using the fundamental theorem of arithmetic, find that the ground state fluctuation vanishes 
t exactly for all excitations. The use of the standard canonical or grand canonical ensembles, on the 

. other hand, gives substantial number fluctuation for the ground state. This difference between the 

microcanonical and canonical results cannot be accounted for within the framework of equilibrium 
00 ' statistical mechanics. 

<n : 
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^ \ After the experimental discovery of BEC in a trapped dilute gas at ultra-low temperatures, much attention has 

S been paid to the problem of number fluctuation in the ground state of the ideal system |2j ||, ||, |[ ^, 0, |[ ^] , as 
well as a weakly interacting bose gas |to[ [H], [l2[ |lj, |l5|, [ll| [l7|. There are also a few papers on ground state 
fluctuations in a trapped fermi gas fig] . There are several reasons for this interest. In standard statistical mechanics, 
number fluctuation is related to density-density correlation, and to the compressibility of the system in the grand 
canonical ensemble (GCE) jnj. Light scattering cross section off the medium, in principle, may be related to the 
ground state fluctuation in the system |p9|| . The so called grand canonical catastrophe in an ideal bose gas, where 
the fluctuation diverges at low temperatures, was already known |2(J. Therefore a more accurate treatment of the 
I | problem was needed for trapped gases. In the microcanonical treatment of number fluctuation from the ground state 
' O ■ in a harmonic trap, the problem is closely related to the combinatorics of partitioning an integer, and thus there was 
an interesting link to number theory It turned out that the result for the ground state number fluctuation was 
very sensitive to the kind of asymptotic approximations that are made. Another aspect that drew much attention in 
the literature was the difference in the calculated results for fluctuation using the canonical and the microcanonical 
formulations (?], |j| . It was pointed out by Navez et al. that for a trapped bose gas below the critical temperature, 
the microcanonical result for fluctuation could be obtained solely using the canonically calculated quantities, which 
■s^j- . in turn may be obtained from the so called Maxwell Demon ensemble [^J, to be explained later. 

In this paper, we give an example of a quantum spectrum that has no number fluctuation in the ground state for 
any excitation energy in the microcanonical ensemble, as a direct corollary of the fundamental theorem of arithmetic. 
The canonical ensemble, on the other hand, yields a dramatically different ground state number fluctuation. The 
method of Navez et al. fails to account for the difference between the microcanonical and canonical results in this 
example. This failure of the canonical (grand canonical) ensembles is due to the peculiar nature of the single-particle 
1""-; . spectrum in our example. Generally, when a large excitation energy is supplied to a system, there are a very large 
number of distinct microscopic configurations that are accessible to it. All these different microstates describe the 
same macro-state of a given excitation energy. The classic example is that of bosons in a harmonic trap, where 
the number of partitions of an integer number, corresponding to the number of microstates, increases exponentially. 
We use, on the other hand, another example from number theory, to propose a system where the excitation energy, 
no matter how large, is locked in one microstate. Consequently, although it is possible to explicitly calculate the 
. canonical or grand canonical partition functions and therefore the thermodynamic entropy for this example, it does 
not approach or equal the information theoretical entropy that can be exactly calculated using number theory [p2[ . 
So far as we know, this constitutes a novel and new example that links number theory and thermodynamics, and 
gives radically different results for the microcanonical and canonical ensembles. 

We consider bosons in a hypothetical trap with a single-particle spectrum (not including the ground state, which 
is at zero energy) 

e. p = hip , (1) 

where p runs over the prime numbers 2,3,5, ... Of course, such a spectrum is not realizable experimentally, and it 
is merely a means of performing a thought experiment. We shall use, in what follows, both a truncated sequence of 
primes, as well as the infinite sequence when we perform the canonical calculations for fluctuation. First, however, we 
perform the exact calculation for number fluctuation from the ground state. Suppose that there are N bosons in the 
ground state at zero energy, and an excitation energy E x is given to the system. In how many ways can this energy 
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be shared amongst the bosons by this spectrum ? Before giving the answer, we remind the reader of the fundamental 
theorem of arithmetic, which states that every positive integer n can be written in only one way as a product of prime 
numbers |2j| : 

n = P Tp¥-P n r r -, (2) 

where p r 's are distinct prime numbers, and n r 's are positive integers including zero, and need not be distinct. It 
immediately follows from Eq. (||) that if the excitation energy E x — Inn, where the integer n > 2, there is only one 
unique way of exciting the particles from the ground state. If E x ^ Inn, the energy is not absorbed by the quantum 
system. Since the number of bosons excited from the ground state, for a given E x , is unique for this system, the 
number fluctuation in the ground state is identically zero ! Moreover, this conclusion is valid whether we take in 
Eq. (Q) an upper cut-off in the prime number, or the infinite sequence of all the primes. The information-theoretic 
entropy at an excitation energy E x is 

S{E X ) = -^T^mP, , (3) 

i 

where Pi is the probability of excitation of the microstate i. Since only one microstate contributes with unit probability, 
and all others have Pi = (when E x = Inn), the entropy 5 = 0. It is also straight forward to calculate the ground 
state population iVo as a function of the excitation energy E x . For this purpose, we truncate the spectrum given by 
Eq. (Q) to the first million primes, with a cutoff denoted by p*, and take N = 100. We shall display the numerical 
results after describing the canonical calculations. 

Our next task is to calculate these quantities in the canonical ensemble, and see if the differences in the mi- 
crocanonical and canonical results may be accounted for (as N — > oo) using the method of Navez et al. We 
first do the calculation for the truncated spectrum. The one-body canonical partition function is then given by 
Z\(0) — 1 + Sp=2 ex P( — /3mp). The N-body bosonic canonical partition function is obtained by using the recursion 
relation |2l 



Zn ^ = Jf E Zi(*0)Zn-,(0) • (4) 

s— 1 

Once Zff is found, the ground state occupation N — {n ) N and the ground state number fluctuation for the canonical 
ensemble can be readily be com puted HI ||. We define 5N$ = {(n$) N ~ (n Q ) 2 N ), where the RHS is calculated 
using Eqs. (24) and (25) of Ref.([18)). In Fig. la) and b), we display the results of the canonical calculations for the 
ground state occupancy fraction Nq/N and the ground state fluctuation (SNq) 1 ^ 2 /N for N = 100 as a function of 
temperature T with the truncated spectrum of the first million primes. For comparison, we also show the results of 
the corresponding grand canonical calculation. The grand canonical catastrophe for the number fluctuation is clearly 
evident. It is also easy to calculate the canonical (equilibrium) entropy S = ]xiZpf(/3) + (5E X . The comparison with 
the combinatorial (or microcanonical) results requires that we convert the canonical temperature to excitation energy 
E x . This is done by using the standard relation E x = ~ 91n ga f ■ In Fig. 2, the canonical fractional occupancy of the 
bosons in the excited states, (N e ) /N, for N — 100, is compared with the (exact) combinatorial (or microcanonical) 
calculation as a function of the excitation energy E x . Although the canonical and the grandcanonical (N e ) /N are 
nearly identical, the corresponding microcanonical quantity is radically different. This anomaly persists even as 
N — » 00, showing the breakdown of the equivalence between the microcanonical and the other ensembles. In Fig. 3, 
we display the behaviour of the canonical entropy as a function of T and E x . The microcanonical entropy S(E X ), of 
course, is zero, and so is the number fluctuation 5N{j. Thus we find that the canonical results have no resemblance 
with the exact microcanonical ones. All these calculations were performed for a truncated spectrum (first million 
primes) of lnp, as specified earlier, and for TV = 100. 

We shall now use the procedure of Navez et al. |Q to check if the microcanonical results may be obtained from a 
canonical calculation as N — > 00. These authors constructed the so called Maxwell's Demon ensemble in which the 
ground state (for T < T c ) was taken to be the reservoir of bosons that could exchange particles with the rest of the 
subsystem (of the excited spectrum) without exchanging energy. Denoting the grand canonical partition function of 
the excited subsystem by 5 e (a,/3), with a = /3fi, it was shown that the canonical occupancy of the excited states, 
(N e ), and the number fluctuation (S 2 N e ) could be obtained from the first and the second derivative of S e with respect 
to a, and then putting a = 0. It was further noted that the microcanonical number fluctuation for the excited 
particles was related to the canonical quantities by the relation 

/x2 N \oo _/,2 N] oo [( SN eSE)^ N } 2 

\° JV e/MC - \° Me/CN (S 2 E)^ ' ^ ' 
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where the superscript oo denotes N — > oo. This worked beautifully for harmonic traps in various dimensions. These 
calculations, for our system, are also easily done for (3 > 1. We now consider the spectrum (]]]) to be the infinite 
sequence of the primes, and evaluate the RHS of Eq. (^). We readily obtain the convergent expression (for j3 > 1) 

(lnp)p' 3 ] 2 

(S N e )f IC = - — {lnp)2pl3 - ■ (6) 

p yF ' t-^p (pP-i) 2 

The RHS of Eq. (0) is nonzero, and therefore does not agree with the microcanonical result. The failure of the 
above formalism of Navez et al. M is because of the very special nature of the single-particle spectrum ([!]) . Consider 
constructing the N-particle canonical partition function Zn(0) from this spectrum, as we had done. As N — > oo, a 
little thought will show that — ► £(/?), where C(/3) = ^ n -\ is the Riemann zeta function. This is because we 
are allowed to span over all E in calculating Zn((3). Similarly, the grand partition function S e (a,/3) with a = is 
none other than the Euler product representation of the Riemann zeta function |2^] . This has a density of states 
growing exponentially with E, and has been studied in connection with the limiting hadronic temperature p6| . By 
contrast, the number of accessible states with a fixed N and E in the isolated microcanonical set-up does not increase 
at all. Since the energy remains locked in one microstate, the isolated system cannot be described through the usual 
concepts of statistical mechanics. 

This research is supported by NERSC (Canada). The authors would like to thank David M. Cooke and Jamal 
Sakhr for helpful discussions. 
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FIG. 1: a) Average occupancy in the ground state No/N versus temperature T for N = 100 in the canonical and grand 
canonical ensembles, b) Plot of the relative ground state number fluctuation in both ensembles. Note the steep rise in the 
grand canonical fluctuation. 



FIG. 2: a) Plot of the average occupancy in the excited states N ex /N , for N = 100, versus the excitation energy E x in the 
canonical ensemble (continuous bold curve), compared with the exact microcanonical calculation. The latter calculation is done 
for E x = In n, where n is an integer, and the results are shown by dark points. These are joined by dotted lines to emphasize 
their zigzag character. For example, the sixth point (including 0) corresponds to E x — In 6, and gives N e = 2, corresponding 
to the prime factor decomposition 2x3. 



FIG. 3: Plot of the canonical entropy as a function of temperature T on the left, and excitation energy E x on the right, for 
N = 100. The microcanonical entropy is zero. 
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